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Abstract 

We compute the hierarchical (^"^-trajectory in terms of perturbation 
theory in a running coupling. In the three dimensional case we resolve 
a singularity due to resonance of power counting factors in terms of 
logarithms of the running coupling. Numerical data is presented and 
the limits of validity explored. We also compute moving eigenvalues 
and eigenvectors on the trajectory as well as their fusion rules. 



1 Introduction 



In the block spin renormalization scheme of Wilson [Wil71, KW74] renormal- 
ized theories come as renormalized trajectories of effective actions. Departing 
from a bare action the renormalized trajectory is reached by an infinite iter- 
ation of block spin transformations. For this limit to exist the bare couplings 
have to be tuned as the number of block spin transformations is increased. 
Consider an asymptotically free model at weak coupling. There the point is 
to keep couplings under control which increase by value under a block spin 
transformation. Such couplings are called relevant. In weakly coupled models 
they can be identified by naive power counting. This renormalization scheme 
has been beautifully implemented both within and beyond perturbation the- 
ory. Let us mention the work of Polchinski [Pol84] , Gawedzki and Kupiainen 
[GK84] , Gallavotti [Gal85], and Rivasseau [Riv9I] as a guide to the extensive 
literature. The underlying picture of an ultraviolet asymptotically free model 
is to think of the renormalized trajectory as unstable manifold of a trivial 
fixed point. Although this picture has been in mind behind block spin renor- 
malization since the very beginning [KW74] it has not been formalized yet 
to an approach free of a bare action. This paper is a contribution to fill this 
gap. It extends the analysis begun in [WX94] and [WX95] in the context 
of renormalization group improved actions for the two dimensional 0{N)- 
invariant nonlinear cr-model. Here we will work it out for the ^"^-trajectory 
in the hierarchical approximation. The hierarchical model was invented by 
Dyson [Dys69] and Baker [Bak72] and has enjoyed the attention of Bleher 
and Sinai [BS73], Collet and Eckmann [CE78], Koch and Wittwer [KW9I], 
Felder [Fel87], and Pordt [Por93], to mention a few. The ^"^-trajectory will 
be defined as a curve which departs the trivial fixed point in the ^"^-direction. 
Technically we perform a renormalized perturbation expansion in a running 
coupling. In the three dimensional case we perform a perturbation expan- 
sion in a running coupling and its logarithm. The dynamical principle which 
proves to be strong enough to determine the trajectory at least in perturba- 
tion theory is stability under the renormalization group. With stability we 
mean that the trajectory is left invariant under a transformation as a set in 
theory space. Recall that a renormalized action always comes together with 
a sequence of descendents generated by further block spin transformations. 
Fven in the case of a discrete transformation this sequence will prove to con- 
sist of points on a continuous curve in theory space which is stable under the 
renormalization group. It is the computation of this curve we address. The 
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result is an iterative form of renormalized perturbation theory. Its closest 
relatives in the literature are the renormalized tree expansions of Gallavotti 
and collaborators [Gal85, GN85a, GN85b]. A pedagocial exposure can be 
found in [FHRW88] . Our expansion is however free of divergencies piled up 
in standard perturbation theory by infinitely iterated renormalization group 
transformations from the very beginning. Surprisingly we do not need to 
treat relevant and irrelevant couplings on a different footing. It will involve 
neither bare couplings nor renormalization conditions in the original sense. A 
renormalization group transformation in our approach translates to a trans- 
formation of the running coupling according to some /3-function. We will 
consider in particular a choice of coordinate whose associated /3-function is 
exactly linear. This idea has also appeared in [EW84] and references therein. 
Renormalized perturbation theory furthermore will surprise us with a se- 
quence of discrete poles at special rational dimensions. These poles will be 
traced back to certain resonance conditions on the scaling dimensions of pow- 
ers of fields. In particular the case of three dimensions will be shown to be 
resonant. We will resolve the associated singularity by a double expansion 
in both the running coupling and its logarithm. The expansion will then be 
extended to the computation of moving eigenvectors in the sense of [WX95] 
on the renormalized trajectory and their fusion rules. Finally we perform a 
numerical test of our renormalized actions. As expected they prove to work 
well in a small field region. The extension of our program to full models 
is under way. A prototype with momentum space regularization has been 
developed in [Wie]. 

2 Hierarchical renormalization group 

The hierarchical renormalization group in the form advocated by Gawedzki 
and Kupiainen [GK84] is a theory of the non linear transformation 





on some space of Boltzmann factors Z[(f>). In the scalar theory ^ is a single 
real field variable. 
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is the Gaussian measure on IR with mean zero and covariance 7. The pa- 
rameters of (1) are the Euchdean dimension D and the block scale L. The 
subspace of even Boltzmann factors Z[ — (f>) = Z[(f>) is stable under (1). We 
will restrict our attention to this subspace. Let the potential be given by 
Z[(f>) = exp( — The transformation for the potential is 

nv{i^) = -L^\og[j dfi.io exp(-y(zi-fv + C))) • (3) 

The below analysis will be done in terms of the potential. The method will be 
perturbation theory. The question of stability bounds will not be addressed. 
Regarding mathematical aspects of (1) and (3) we refer to the work of Collet 
and Eckmann [CE78], Gawedzki and Kupiainen [GK84] and of Koch and 
Wittwer [KW91]. 



3 The trivial fixed point 

(3) has a trivial fixed point V^{(f>) = 0. This fixed point is the hierarchical 
massless free field. The linearization of (3) at this trivial fixed point is given 
by 



V' + C • (4) 



This linearization is diagonalizable. The eigenvectors are normal ordered 
products 



•7/ ■ 



df 



^_^exp(j> —] (5) 



with normal ordering covariance 7/ = (l — L^~^^ ^ 7. The normal ordering 
covariance has been chosen in order to be invariant with respect to integration 
with dfi^. Its singularity at D = 2 is an infrared singularity of the hierarchical 
massless free field in two dimensions. The eigenvalues are 

A„ = Z^+Ki-f ) . (6) 

The eigenvalue of : c/)'^ is A4 = L^~^ . The eigenvector : cj)^ is therefore 
relevant for D < 4, marginal for D = 4, and irrelevant for D > 4 dimen- 
sions. Perturbation theory can be used to compute corrections to (4) in a 
neighbourhood of V^{(f>). 
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4 The (j)^ - trajectory 

Let us define a curve V{(f>^g) in tlie space of potentials parametrized by a 
local coordinate g. We call it the (p'^ - trajetory. We expand the potential 

oo 

V{<l>,g) = Y,V2n{9) (7) 

n=0 

in the base of eigenvectors (5). A natural coordinate in a vicinity of V^{(f>) 
is the (p'^ - coupling defined by V4{g) = g. Let us use it for a moment. Let 
the cf)'^ - trajectory then be the curve V{(f>^g) defined by the following two 
conditions: 

1) V{(f>^g) is stable under TZ. Then there exists a function l3[g) such that 

nV{<l>,g) = V{<l>J{g)) . (8) 

The function l3[g) is of course coordinate dependent. With the (p'^ - 
coupling as coordinate it is called /3 - function. 

2) V{(p^g) visits the trivial fixed point V^{(p) at g = 0. The tangent to 
V{<p,g) at V^{<p) is given by 

d 



dg 



V{<l>,g) = :<l>^:,, (9) 

3=0 



This condition is equivalent with V4{g) = g + 0{g^) together with 
V2n{g) = 0{g'), n^2. 

The (p'^ - trajectory is the object of principal interest in massless (p'^ - theory 
at weak coupling. 

5 Perturbation theory 

The (p'^ - trajectory can be computed by perturbation theory in g as solution 
to (8) and (9). Potentials on the (p'^ - trajectory are said to scale. A potential 
V{(p^g) is said to scale to order s in g if there exists a function 

m = ^^'Ha) + 0{g'+') , 

s 

^^'\9) = Y.Kg^ , (10) 

r = l 
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such that 



V{<t>,g) = V^^\<t>,g) + 0{g^^') , 

nV^'\c^,g) = V^'\c^Ji{g)) + 0{g'^') , (11) 

and 

V^'\^^,g)=g:^':„ . (12) 

The scheme is to compute li^^'^^\g) and V^^'^^\<f), g) given fi^'\g) and 
V^^\<f),g) to some order s. Let us explain it in some detail at the case of 
D = 4: dimensions, block scale Z = 2, and covariance 7 = 1. Then the normal 
ordering covariance is 7/ = |. Computing a block spin transformation (3), 
we speak of V{(f>) as bare potential and of TZV{(f>) as effective potential. The 
point of departure is (12). Anticipating the terms generated in TZV^^\(f)^ g) 
to second order in g we make the ansatz 



V^'H^^g) = Cog' + C2g' ■.^'■.+g:^': +ceg' : f : . (13) 

The coefficients are determined by the condition that (11) be fulfilled to 
second order. (13) is mapped to 

^440 

nv('^{<i>.g{g')) = (16co - ^)/ + (4c2 - 448)/ : : 

+9' : : +( J - 2)5" : f : +0(/) . (14) 

Here the effective coupling defined as the coefficient of : : in the effective 
potential is given by 

g'{g) =g- GO^/^ + 0(/) . (15) 

Comparing the effective potential as a function of the effective coupling with 
the bare potential as a function of the bare coupling we conclude that 

1088 448 8 

CO = ^ , c, = — , ce = -- (16) 

on the cf)'^ - trajectory. The coefficients of the /3 - function (10) to this order 



are 



61 = 1 , 62 = -60 . (17) 
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It follows that g is marginally irrelevant in four dimensions. This completes 
the first step. It is iterated in the obvious manner. The general form of the 
order s approximation is 



2n 



s 

r = 2 



ci'Kg) = g 



-2)1 



(g) = ^ C2n,rg' , n ^ 3 . (18) 



r=n — l 



It includes all normal ordered products generated in the effective potential by 
(3) from (12) to order s in g. The iteration proceeds as above with the order 
s + 1 ansatz of the form (18). The condition (11) yields a system of linear 
equations for the order s + 1 coefficients. (To highest order the coefficients 
have no other choice.) This system has a unique solution: the (p'^ - trajectory. 
Note that the coefficient bs+i of the /3 - function is already determined by 
V^'^(l),g). For inst ance (15) does not contain any of the coefficients in (13). 
The expansion can be computed to higher orders using computer algebra. 
To third order we find 

(3), , 1088 2 54784 . 

Co {g) = -^g — ^g , 

cP(5) = -^/ + 352/, 

cfig) = f / (19) 



together with 



/3(3)(5) =5 -60/ + 8880/ . (20) 



Let us remark that the perturbation coefficients (10) and (11) come with 
alternating signs. The coefficients show a frightening increase in absolute 
value with the order in g. The full series is not expected to converge. Note 
that the coefficients look better when g is replaced by g/4:\. 
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6 Resonances 



We can apply the above scheme to compute the (p'^ - trajectory in other than 
four dimensions. The solution is again of the form (18). We do however 
encounter a new phenomenon. The improvement coefficients C2^^ exhibit 
poles at certain discrete points as functions of the dimension parameter D. 
We call this phenomenon resonance because it can be traced back to the 
fact that certain scaling parameters become powers of one another at these 
points. This happens particular in three dimensions. Let us consider the 
transformation (3) with block scale L = 2 and covariance 7 = 1, but this 
time with arbitrary value of D. We express the dependence on D in terms 
of a variable a = . To third order in g the improvement coefficients are 
given by the rational functions 

(3) _ 12a3(a + 4)(a^ + 16) ^ 
" (a3 - 256)(a -4)3 ^ 



288a\a^ + 32a^ + 512a^ + 4096)(a + if 3 



2 

19^ 



[a + 8)(a - 8)(a2 + 64)(a3 - 256)(a - 4) 
(3) _ 48a2(a2 + 4a + 16) 2 
" (a-8)(a + 8)(a-4)2^ 

384a3(7a^ + 46a^ + 288a'^ + 1728a2 - 1024a - 22528) 
(a + 8)(a - 8)(a3 - 1024)(a - 4)3 



,3 



(3) 8 2 576a(a + 6) 3 

" "F " (a_4)(a-64)^ ' 
(3) 32 3 

4 = y5 • (21, 
The /3 - function to this order in g is given by the function 



16 576(a + 4) , 

—9 — r — 

a a[a — 4j 



256(215a4 + 1400a3 - 10128a ^ 
^ a(a - 8)(a + 8)(a - 4)2 ^ 
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256 95744a + 355328 3 , , 

;9 ■ (22) 



a(a - 8)(a + 8)(a - 4) 



2- 
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Singularities appear at positive dimensions 





g' 




a 


4 2563 


8 


1024^ 64 


D 


2 ^ 

^ 3 


3 


^ 6. 

3 



As one goes to higher orders in g more and more poles show up in (21) and 
(22). Let us have a closer look at the three dimensional pole to second order 
in g. Inserting (13) into (3) at three dimensions we find 

nV^'\<P,g{g')) = (2co - 360) + (C2 - 336)/ 1(^2: 

+g' : : + (I - 2) g" : : +0{g") . (24) 

The parameters of (3) are D = 3, Z = 2, and 7 = 1. The normal ordering 
covariance is 7/ = 27. The /3 - function to this order is given by 

g'{g) = 2g- 21Qg' + 0{g^) . (25) 

The : (p'^ : - coupling is therefore relevant in three dimensions. From (24) we 
would conclude that 

Co = 360 , C6 = . (26) 

But there exists no solution to the equation for C2 (besides infinity). It follows 
that the : : - coupling cannot be written in terms of a power series in the 
: (p'^ : - coupling on the (p'^ - trajectory in three dimensions. The point is that 
the : cf)^ : - coupling fiows like 

336 

- 336n22'^ = --^log(5)5^ (27) 
log(2) 

with g = 2"^ upon iteration of (24). This suggests a double expansion in both 
g and log(5'). 



7 Linear l3 - function 

So far we have used the (p'^ - coupling as coordinate for the (p'^ - trajectory. It 
was defined by the condition V4{g) = g following (7). This coordinate leads 
to unnecessary complications when dealing with its logarithm. A better 
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coordinate which is also interesting by itself is the linear coordinate defined 
by the condition that the /3 - function be exactly given by 



(28) 



on the (p'^ - trajectory. In terms of this linear /3 - function the transformation 
on the (p'^ - trajectory looks exactly like the linearized renormalization group. 
The definitions (8) and (9) remain untouched by (28). Using the expansion 
with the linear /3 - function 



also becomes a computable power series. The strategy is the same as above. 
For the case of L = 2 and 7 = 1 the (p'^ - trajectory is given by 



12a3(a + 4)(a2 + 16) ^ 
(a3 - 256)(a -4)3 ^ 

576a4(a3 + 8a2 + 8a + 256)(a + 4)2 3 
^{a- 16)(a + 8)(a - 8)(a2 + 64)(a - 4)^^ 

48a2(a2 + 4a + 16) 2 
(a-8)(a + 8)(a-4)2^ 

768a3(a6 + 69a^ + 368a* - 2880a3) 
^{a- 16)(a + 8)(a - 8){a^ - 1024)(a - 4)3 

768a3(-22528a2 - 144384a - 475136) 
^{a- 16)(a + 8)(a - 8)(a3 - 1024)(a - 4)3 
36a(a + 4) ^ 
(a-4)(a- 16)^ 
16a2(53a^ - 3336a4 - 24752a3) 3 
~ (a + 8)(a - 8)(a + 16)(a - 16)2(a - 4)2^ 
16a2(149248a2 + 1342464a + 5685248) 3 
~ (a + 8)(a - 8)(a + 16)(a - 16)2(a - 4)2^ 

8 2 384a(2a2 - 45a - 272) 3 
~3^ ~ (a - 4)(a - 16)(a - 64)^ ' 



(29) 



r = 2 




32 



(30) 
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To each order of perturbation theory we find a system of hnear equations 
which has a unique solution. The coefficients are again rational functions in 
a = with poles at resonant dimensions. Note that (30) has an additional 
pole at D = 4 as compared with (21). In four dimensions (28) becomes the 
identity and (8) a fixed point equation. Resonances are now easily under- 
stood. To order m in g the transformation (3) acts as 

^^D+n(2-D)-m(A-D)^^^g^m , ^2n . _ ^3^^ 

A resonance thus occurs if 

D + n(2 - D) - m(4 - D) = . (32) 

In the case of D = 3 dimensions this condition becomes 

3 - n - m = . (33) 

Since m ^ 2 in this business we find only two resonant terms (n, m) in three 
dimensions: (1,2) and (0,3). The former is a mass resonance, the latter a 
vacuum resonance. The resonant dimensions are rational and given by 

„ 4m — 2n , ^ 

D = . 34 

1 — n + m 

In particular table (23) is immediatly reproduced. An interesting variation 
of the linear /3 - function consists of replacing (28) by 

f^is) = big + b2g^ , 
6i = Z^-^ 



7 

2 I tD 



b, = 36Z-^|^, (35) 

truncating the /3 -function (8) to second (or more generally to any fixed) order 
of perturbation theory in the (p'^ - coupling g. This /3 - function has a fixed 
point at finite g for D < 4. Supposing that our expansion was summable 
in some sense it follows that the nontrivial fixed point of Koch and Wittwer 
is on the (p'^ - trajectory. Recall that in these coordinates (8) holds exactly 
with /3 - function (35). It follows that a fixed point of (35) is a fixed point of 
(3). We will not pursue this line of thought further at this instant. The /3 - 
function (28) is ideally suited for the double expansion in both g and log[g). 
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8 Perturbation theory in g and log[g 

Let us consider the (p'^ - trajectory in D = 3 dimensions, again with L = 2 
and 7 = 1. The problem that C2 cannot be determined in (24) such that it 
remains invariant to second order can be cured as follows. We use the linear 
/3 - function defined by (28) and expand the potential in both g and 

« = log(5) • (36) 

In the expansion we treat k, as an independent variable having the same order 
as Note however that our expansion consists of at least differentiable 
combinations of g and k. To second order in g^ we can replace (13) by the 
ansatz 

y(2)((^, (/) = Cog' + (C2 + C2,i/€)(/2 : +{g + eg') +ceg' : f : . (37) 

Here we have left out all terms which anyway turn out to be zero on the (f)'^ - 
trajectory. In terms of g' = 2g and k' = /€ + log(2) the ansatz (37) is mapped 
by (3) to 

nV^'H.^,9{/)) = (2CO- 360)/ 

+ (C2 - log(2)c2,i - 336 + C2,i/€') g'' : <!>' : 

+ [g' + (I - 54) g") : : 

+ (I - 2) g" : : +0{g") . (38) 

As a consequence (37) reproduces its form up to a change of the running 
coupling (28) if and only if 

8 336 
Co = 360 , C4 = -108 , C6 = -- , C2,i = -- — — . (39) 

3 log(2) 

The parameter C2 is free. To second order in g we find a one parameter family 
of solutions to (8) and (9). The free parameter is associated with the mass 
resonance (1,2) of (33). One immediately anticipates another free parameter 
to third order in g coming with the vacuum resonance (0,3). This is indeed 
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the case. The general solution of (8) and (9) to third order in g is given by 

y(3)(<^,/ = 360/ + ^A + co/ 

log(2) 

/ , 336 , , , . 12096 . 

+ - + 116928 - 36c2 / + 

V log(2) log(2) 

f 9 / 8 , . 896 . 

+ (g- 108(/2 + 17520 - -C2 / + 

V 3 log(2) 

+ (-^^ + 864/) 

+ f / : f : (40) 

including Cq and C2 as free parameters. Thereafter we have no further free 
parameters in the higher order coefficients. Thus we have to supplement 
the definition of the cf)'^ - trajectory given by (8) and (9) by two additional 
conditions on Cq and C2 in three dimensions to single out a curve in the space 
of interactions. We choose 

Co = , C2 = . (41) 

With this choice the perturbation theory has a minimal number of vertices. 
(41) can be thought of as additional renormalization conditions. 
At higher orders the scheme explained above iterates. The general form of 
the potential at order s is 

V^'H<P,g) = Y.^t\9) ■■ <P'^ ■■ , 

n=0 

r = 2 t=0 

s [r/2] 

r = 2 t=0 

yt\9) = E 5]^2.,.,t/«* , n ^ 3 . (42) 

r=n-l t=0 

Here the third order coefficients are given by (40) and , for instance, (41). To 
each further order of perturbation theory we meet a system of linear equations 
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possessing a unique solution lor the coefficients in (42): the (p'^ - trajectory 
in terms ol the double expansion. Recall that k, should be substituted by 
log{g) in (42). 



9 Numerical calculation of the renormalized 
trajectory 

Hierarchical renormalization group flows can also be computed using stan- 
dard numerical methods. Let us perform a numerical investigation of the 
transformation (3) in order to determine the limits of validity of the expan- 
sion (42) in g and log[g). Furthermore we want to investigate the large field 
behaviour of our expansion. We choose the following numerical setup. To 
calculate the transformation (3) iteratively we sample the potential V at N 
equidistant points betwen and cf>max- Then we perform a cubic spline in- 
terpolation and integrate using standard NAG library functions. To reduce 
the error due to boundary effects at ^ = cf>max we always choose cf>max so that 
V{(i),max) = V,max with V,max large cuough; for example V,max = 20 will do. For 
^ > <t^max we set 



and choose a and h so that the first derivative of V at cf>max is continuous. 
Here with Vht we mean the quadratic high temperature fixed point of (3). 
Field asymptotics have been investigated by Koch and Wittwer in their work 
on the nontrivial double well fixed point [KW91]. That means that we sup- 
plement our numerical calculation with the expectation that the potentials 
on the renormalized trajectory have HT-like asymptotic behaviour. We let 
the fiuctuation field C, vary between —C,max and C,max and choose C,max = 20. 
This produces errors which can be neglected. As in the case of the expan- 
sion we restrict our attention to the room of symmetric potentials which is 
invariant under (3). All potentials were calculated with iv = 2, _D = 3, 7 = 1 
and = 401. To determine the renormalized trajectory according to our 
definition in section 4 supplemented by the condition C2 = 0, see section 8, 




(43) 
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we proceed as follows. We start with a bare potential 



V 




4 



(44) 



and choose go to be a sufficiently small number, that means that an iteration 
starting with L~'"^ go would yield the same trajectory if m G N. We have 
chosen go = 10~^. Note that we have to choose the third order of our 
perturbation expansion as input because the second order wouldn't be stable 
under iterated renormalization group transformations because of the large 
field behaviour. It is convenient to normalize the potential with the condition 
V{cf> = 0) = 0. The numerically determined potential will be refered to as 
exact in the following. We perform 12 iterations of the transformation (3). 
The analogous perturbative potentials can be found by solving the equation 



for §0 (which gives approximately 10 ^ of course) and then follow the renor- 
malization group fiow to 



For the calculation of the perturbative potentials we used the seventh or- 
der of our perturbation expansion. The comparison between the exact and 
the perturbative potentials can be seen in the following pictures. The ex- 
act potentials correspond to the continuous lines the perturbative potentials 
correspond to the dashed lines. In figure (1) one can see that after 8 renor- 
malization group steps the exact and the perturbative data are nearly the 
same. If however the number of renormalization group steps exceeds 9, there 
is a clear digression for cf) ^ 18. This can be seen in figure (2). That means 
that the large field behaviour of the perturbative expansion is wrong, since 
the asymptotic behaviour of the exact potential is quadratic. If the number 
of renormalization group steps is bigger than 10 the perturbative expansion 




(45) 



L^go , n = 0, . . . , 13 . 



(46) 
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is only valid for very small fields. We therefore conclude that our expansion 
is not valid any more for g > 10""^ . The borderline has been illustrated in 
figure (3). One can however improve the perturbative data if one is willed 
to work with pade approximants. In figure (4) we compare the (9, 7) pade 
approximant of our perturbation series in the variable cf) ( which has the 
expected asymptotic behaviour ) with the numerical data after 11 renormal- 
ization group steps. Amazingly enough both curves coincide. We will not go 
in details here. Let us only mention that some pade approximants develop 
additional unphysical poles at real values of cf). 



10 Observables 

Consider a local observable 0{cf>{x)) in a full theory . In the hierarchical 
approximation it corresponds to an observable function 0{cf>),cf> G IR. Ob- 
servables are transformed according to the linearized renormalization group. 
The linearization CyTl-O of the transformation (3) in the direction O at the 
potential defined by 



£y7^(9(V') = |- 

OZ 



n{V + zO){i^) (47) 

z=0 



is given by 

r ^n(i^ / dniC) o{L'-T^p + C) exp(-y(X^-f v> + 0) 

LvrCUitp) = ^ ^ 

/J/.,(C) exp(-y(Zi--V + C)) (48) 

after division by . Note that this is analogous to the linear block spin 
transformation of local observables in the full model. Now we want to find 
the running eigenvectors and eigenvalues of CyTl-O along the renormalized 
trajectory. Here an observable is called a running eigenvector if it satisfies 
the equation 

Cv^,nO{<t>, g) = e{m)0{<t>. m) ■ (49) 

e[g) is the corresponding running eigenvalue. Note that O and e are para- 
metrized by the perturbative expansion parameter g. /3{g) means the /3 - 
function, for example (28). We already know the perturbative expansion of 
the potential V along the renormalized trajectory (in a certain neighbour- 
hood of the trivial fixed point) so all we have to do is to solve (49) for O 
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perturbatively order by order. We start our iteration with the eigenvectors 
: at the trivial fixed point, see section 3. The improved observable to 

the order s of perturbation theory which equals : -.^i for g = will be 
called On \(/),g). The corresponding running eigenvalue will be denoted by 
en \g). Let us explain the scheme of the iteration to some detail by calculat- 
ing (92^^(^, (/) in three dimensions, using the linear /3-function. To first order 
of perturbation theory for this observable we expect 

^2 \(t^^9) = CAflg : (j)'^ : + : (j? : +(co,o + CQ^iK)g . (50) 

We now calculate an effective observable following the transformation (48) 
to get 

Cv,,nO^^\^^,g) = Qc4,o - l) 5 : : 

+ Qco,o + ^Co,i {k - ln(2))^ g . (51) 

Therefore we conclude that the running eigenvalue £^2 ^9) ^^^^ order is 
given by 



e 



2 



(5) = t: - 95 . (52) 



Now we divide by £^2 ^9) ^^'^ 



2 



Oi:ijj = |^-C4,o - 2J , : 
+ :^': 

+ (co,o + co,i(/€-ln(2)))5 , (53) 

after scaling of g. That means that we have the desired invariance if we set 

g 

C4,o = ~3 ' "^o,! = . (54) 

Analogous to the determination of the potential along the renormalized tra- 
jectory we find that one parameter , Co,o, can be freely chosen. This can be 



16 



viewed as an extra renormalization condition which is imposed on the observ- 
able 02. In the general case this occurs under the following circumstances: 
Let us have a look at a term eg' : i/j"^ : of the observable On to l-th order of 
perturbation theory. This is mapped to 

cZ™(i-f )+'(^-4)-Ki-f : ^™ : (55) 

because we divide by e„((/). That means that resonances occur if 

D = 56 



or 

2/(4-1)) 



m 



D 



+ n. (57) 



For the case D = 3 this means that we have extra renormalization conditions 
whenever 

m = n - 2/ . (58) 

So for the observable O4 for example we have two such conditions. 

In the following table we display three observables using the parameters 

iv = 2, _D = 3,7 = 1. Here and in the rest of the article we choose all free 
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parameters to be zero. 
O^i^.g) = 1, 



7/ 

4 , -, . . 2 or.. 



+ ( -35 + 14%' - 896^ - 23360/ ^ • 
+ ( —r- - 1728/ / • r -7' 



•7/ 



,3 



256/ ^7 



.^f3i/ , X 9 2016 

0^i'{^,g) = -960 / + ——f^ - 632448 5 



3 



ln(2) 



+ : :7/ 



85 , .00 ,2 1792 //€ 483072/ . 4 



+ ( -^ + 432(/ 
+ 



3 ^ ln(2) 
112/ 164160/ 



•7/ 



640/ 



9 31 

8 . 



9 -v^ 



(59) 



The corresponding eigenvalues are 

(3) V2 84^2/^/ ^ 2 3528 y2/c/ /- 3 

= - + n^ + 28V2, - 1,(2) - 34608^2, , 

= 1/2 - 9, + + 1232/ -IIJ^- 287736/. (60) 

We have calculated the running eigenvectors and the running eigenvalues to 
seventh order. This is used for the fusion rules which will be derived in the 
next section and for the comparison with the numerical data . 



11 Fusion rules 

With the help of the observables On on the renormalized trajectory one can 
compute more general correlation functions. This task requires the introduc- 
tion of a new concept, the hierarchical fusion rules. Given two observables 
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On{(t>',g) and Ojnid^-tQ) these are defined by 

CO 

0^{c^,g)0^{c^,g) = Y,Nl^{g)Oi{c^,g) . (61) 
/=o 

Note that with the coefficient we get a symmetric bihnear form <, > on 
the space of observables, defined by 

< O.M,g),0.r,{^,g) >:=Kjg) . (62) 

In the thermodynamic hmit only the overlap of an observable with the con- 
stant term will survive. 

According to our general concept we expand the fusion coefficients Nl^{g) 
perturbatively: 

s 

Kl:\9)=Y.^n-r.ik)9' ■ (63) 
k=0 

Now it is straightforward to calculate the fusion coefficients to some order of 
perturbation theory. To zeroth order we recover the well known fusion rules 
for normal ordered products ,of course, for example 

(9!°)(<^,5)0!°)((^,5) = 2(9r(<^,5) + 0f(<^,5) , 
(9f(<^,5)0f((^,5) = 8(9r(<^,5) + 8(9f (<^,5) + 0f (<^,5) . (64) 
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Here the parameters are D = 3, 7 = 1, and therefore the normal ordering 
covariance 7/ equals 2. To second order we get for example 



+ i^l-32g + g' (224^^ + 5504) j 0^^\^, g) 

-9Qg'0^i\<P,g) , 
0^^\<P,g)0^^\<P,g) = f8 - 7205 



Note that we have chosen the free parameters to be zero in these formulas. 

12 Numerical calculation of the observables 
and eigenvalues 

Similar to the potentials on the renormalized trajectory we are now going to 
calculate the eigenvalues and observables of our theory in three dimensions 
numerically. The aim is again the determination of the region of validity 



+ {l-16g + 1856g')0^^\^,g) 
-i8g'Of\<P,g) , 
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of our peturbative expansion in g and \og{g) now for the eigenvalues and 
observables. To this end we restrict the action of the linearized renormaliza- 
tion group transformation (48) to the finite dimensional space of observables 
which is spanned by 

, ^ m ^ M . (66) 

After performing an expansion (which is most conveniently done by differen- 
tiation) we get a finite dimensional representation matrix L of the linearized 
renormalization group transformation 

M 

Cv^,n<t>' =Y,U,<t>^ . (67) 

Now we are able to calculate the eigenvalues and observables of L. Because 
the matrix L is of course only an approximation of the transformation (48) we 
have to choose M big enough to get correct results. We have chosen Af = 8 
and expect that the first four eigenvalues and eigenvectors do not suffer from 
big errors due to the trucation. In figure (5) we have plotted the four largest 
eigenvalues of transformation (48) against the number of renormalization 
group steps. The crosses correspond to the exact (i.e. numerical) calculation 
whereas the boxes correspond to our perturbation expansion. We start at 
the perturbative potential at go = 10~^, that means near the trivial fixed 
point. Then we follow the renormalized trajectory performing numerical and 
perturbative renormalization group steps. After six renormalization group 
steps, i.e. at = 2^go perturbatively, compare section 7, one can see the 
first small deviations between the perturbative and the exact eigenvalues. 
After 10 renormalization group steps (i.e. at = 1.0 • 10""^) there is a clear 
distinction between both. We recover and sharpen the former result that 
for g > 10""^ our expansion cannot be said to be valid any more because of 
nonperturbative effects. 

In figure (6) to figure (9) we have plotted the exact and perturbative 
observables O2 and 03 after four and after six renormalization group steps, 
respectively. Here the perturbative data correspond to the dashed lines. For 
large fields we have deviations for both observables which are in some extent 
due to the truncation of the transformation (48). The small field behaviour 
has been replotted in the small pictures respectively in order to illustrate the 
more significant infiuence of the truncation on the observable 03 as well as 
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to illustrate the right behaviour at ^ = of our calculations. Here we just 
recover the approximate normal ordering near the trivial fixed point. Because 
of nonperturbative effects we cannot reach the nontrivial fixed point using the 
naive series for the eigenvalues. It would be very interesting however to find 
a way to sum up our expansion in g and \og{g) in order to calculate critical 
exponents at the nontrivial fixed point. This problem will be investigated in 
our future work. 
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Figure 1: exact and perturbative potential after 8 renormalization group 
steps 
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Figure 2: exact and perturbative potential after 9 renormalization group 
steps 
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Figure 3: exact and perturbative potential after 10 renormalization group 
steps 
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Figure 4: exact potential and pade approximant of the perturbative potential 
after 11 renormalization group steps 
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Figure 5: exact (crosses) and perturbative (boxes) eigenvalues of transfor- 
mation (48) against number of renormalization group steps 
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Figure 6: 02 (^) exact and perturbative after four renormalization group 
steps 
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Figure 8: 03{(f>) exact and perturbative after four renormalization group 
steps 
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